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ABSTRACT 

We formulate the problem of magnetic field dissipation during the accretion 
phase of low-mass star formation, and we carry out the first step of an itera- 
tive solution procedure by assuming that the gas is in free-fall along radial field 
lines. This so-called "kinematic approximation" ignores the back reaction of the 
Lorentz force on the accretion flow. In quasi steady-state, and assuming the 
resistivity coefficient to be spatially uniform, the problem is analytically soluble 
in terms of Legendre's polynomials and confluent hypergeometric functions. The 
dissipation of the magnetic field occurs inside a region of radius inversely propor- 
tional to the mass of the central star (the "Ohm radius"), where the magnetic 
field becomes asymptotically straight and uniform. In our solution, the magnetic 
fiux problem of star formation is avoided because the magnetic fiux dragged in 
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the accreting protostar is always zero. Our results imply that the effective re- 
sistivity of the infalling gas must be higher by several orders of magnitude than 
the microscopic electric resistivity, to avoid conflict with measurements of paleo- 
magnetism in meteorites and with the observed luminosity of regions of low-mass 
star formation. 

Subject headings: ISMxlouds — ISM: magnetic fields — magnet ohydro dynamics 
— planetary systems: protoplanetary disks — stars: formation 



1. Introduction 

If the magnetic field of a collapsing interstellar gas cloud remained frozen in the gas, the 
resulting surface magnetic field of the newborn protostar would exceed observed stellar fields 
by almost four orders of magnitude (the so-called "magnetic flux problem", see e.g., Mestel & 
Spitzer 1956). It follows that the excess magnetic flux of a cloud must be dissipated at some 
stage during the process of star formation. Assuming ideal magnetohydrodynamic (MHD) 
conditions, Galli et al. (2006, hereafter Paper I) obtained an analytical solution for the inner 
regions of an isothermal, magnetized, rotating cloud undergoing gravitational collapse and 
showed that the long-lever arms of the strong field trapped in the central protostar in a 
split-monopole configuration would cause so much magnetic braking as to make impossible 
the formation of a centrifugally supported disk around the central object (see also Allen et 
al. 2003a,b). Thus, the dissipation of magnetic field must occur prior or simultaneosly to 
the formation of a circumstellar disk. In the current paper, we remove the assumption of 
field freezing by allowing the non-ideal effect of finite electric resistivity to operate, and we 
determine the value of the resistivity coefficient required to solve the magnetic flux problem 
during the accretion phase of low-mass star formation. 

The paper is organized as follows: in Section 2 we formulate the basic equations of the 
problem; in Section 3 we estimate the value of the resistivity coefficient in the central region 
of collapse from the available observational constraints; in Section 4 we check the validity of 
our approximations; in Section 5 we compare our results with those of previous works; finally, 
in Section 6, we summarize our conclusions. The reader not interested to the mathematical 
derivation of our results, can go directly to eq. (31). 



2. Formulation of the problem 

Neglecting the effects of ambipolar diffusion (but see Section 6), the evolution of the 
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magnetic field in a collapsing cloud is governed by the induction equation with ohmic dissi- 
pation, 

— + Vx (B X u) = -Vx (r^Vx B), (1) 

where rj is the coefficient of ohmic resistivity, related to the electric conductivity a by 

Assuming axial symmetry, and a purely poloidal magnetic field, we can "uncurl" eq. (1) and 
express it in scalar form as 

— + u.V$ = ^5($), (3) 

where $ is the flux function, defined in spherical coordinates by 



B = Vx 



e 

27rr sin 6 



(4) 



and S is the Stokes operator 



S = — 1 9^ cotg d 

Qj.2 ^2 QQ2 ^2 QQ 



Since B.f. is antisymmetric with respect to the midplane, eq. (4) implies that the flux function 
$ is even with respect to ^ = 7r/2. In addition, because the field is finite on the polar axis 
(except at the origin), the flux function $ must vanish for ^ = and 9 = ti ^ 

Following Paper I, we assume that the velocity is given by free-fall on a point mass M^ 
at the origin. 



u r 



and we ignore the reaction of the magnetic field on the flow (the validity of this assumption 
is to be checked a posteriori). Because the dissipation is likely to occur in a region of space 
across which the flow time is short in comparison with the evolutionary time of the collapse 
of the molecular cloud core, we may look for a quasi-steady state in which the advection of 
magnetic flux by the flow is balanced by ohmic dissipation, 

V^m- (7) 



\ r J dr 



^The physical flux threading a circle of radius r on the midplane is given by the value $(r, 9 = 7r/2). 
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At large r, the flux function must asymptotically approach that of a split monopole, 

lim ^(r,^) = $41- Icos^l), (8) 

i — >oo 

where 

$^ = 27rX:^G^/^M^ (9) 

is the flux trapped by the central protostar, and the parameter Aj, is the mass-to-flux ratio of 
the central split monopole in non-dimensional units (see Paper I). Note that Ay, is a measure 
of the magnetic flux trapped in the star under ideal MHD condition, and does not correspond 
to the actual mass-to-flux ratio measured in young stars (~ 10^-10^ in the same units). In 
Paper I, we obtained 1 ^ A^ < 4 by connecting the analytic inner collapse solution to the 
ideal MHD numerical models of Allen et al. (2003a,b). In the following, we will use A^ ~ 2 
as a flducial value. 



2.1. Nondimensional variables 

To be speciflc and for simplicity, we assume that 77 is spatially constant and we deflne 
the nondimensional variables 

r = rohmX, $(r, 6) = $^0(x, /i), (10) 

where n = cos 6 and 

is the Ohm radius. We may then justify the assumption of a quasi-steady state if the 

fractional variations of the parameters of the problem, M^ and $^, are negligible over a 

diffusion time 

2 s 

^ohm = ^ - ja^- (12) 

Because advection is balanced against diffusion, the expression (12) is also the time it takes 
to cross rohm at the local free-fall velocity 

w(rohm) = = — — , (13) 

\ ^Ohm / V 

which explains why rohm is deflned as in eq. (11). If we adopt the values rj ^ 2 x lO^*' cm^ 
s^^ and M^ ^ 1 Mq (see Sect. 3), we get tohm ~ 3 yr, which is very short compared to the 
timescale 

tacc = ^ ~ 10' y^' (14) 
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over which we may expect M^ or $^ to have significant variations. Thus, the quasi-steady 
approximation is hkely to be a good one. 



With these definitions, eq. (7) becomes 



dx^ dx dij? 



Setting (f){x,fi) = F{x)G{fi), we reduce eq. (15) by separation of variables to the couple of 
ordinary differential equations 

{1 - fi^)G" + AG = 0, (16) 

where A is a separation constant, and 



^F" + x^/^F'-AF = 0. (17) 



X' 

Differentiating eq. (16) with respect to /x, and setting g{fj,) = G'{fj.), we see that g{fj.) 
satisfies Legendre's equation, 

{l-fi^)g"-2fig' + Ag = 0. (18) 

Solutions of Legendre's equation regular at ;U = ±1 require A = /(/ + 1), with / = 0, 1, 2, . . ., 
and are given by Legendre's polynomials Pi{fi) of order /, i.e. g{fi) oc Pii^jj). Using the 
definition of g and eq. (16), we then obtain G{^) oc (1 — fi^)Pl{fi). Therefore, any solution 
of eq. (16) regular on the polar axis can be written as a linear combination of polynomials 
of degree n = / + 1, 

dP 

GM = C„(l - /i')-f^ = nC„[/xP„_i(/i) - PM], (19) 

d/i 



where C„ are arbitrary constants. Defining 



Gn = 



- 2n-l 1 '/' 
_2n(n-l)J ' 

the polynomials Gn{fi) form an orthonormal set with weight (1 — /i^)~^. 



(20) 



-1 d/i 

GnifJ')G„i{fl) = 6nm, (21) 

L i — /i 

as can be easily shown by a simple integration by parts and using the normalization condition 
of Legendre's polynomials. Since the fiux function is an even function of /i over the interval 
— 1 < yW < 1, the index n takes only even values n = 2,4, . . .. Figure 1 shows the functions 
GniO) for n = 2-10. 
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The boundary condition (8) must be expanded in terms of Gn polynomials, 

oo 

lim 0(x,/i) = 1 - |/i| = y^ fnGnifJ'), (22) 

.T — >no < * 



n=2,4,. 



where /„ are the spectral coefficients. Multiplying both sides of this equation by G'm(;u)/(1 — 
/i^) and integrating from /i = — lto/i = lwe easily obtain, using the orthogonality property 
(21), and integrating by parts, 

/„ = 2C„ /' P^^^{^,) d/x = (-l)"/2-i^i!i^lll!L^C„, (23) 

Jo n{n- l){n-2)U 

where the latter equality follows from a formula in Byerly (1959). 

We now turn to eq. (17), with A = n{n — 1). For each n, this equation has to be solved 
under the boundary conditions 

F„ = at a; = 0, lim F„(a;) = /„, (24) 

X— »oo 

where the constants /„ are given by eq. (23). In terms of the new variables z and i/„ defined 
by 

z = 2y^ and F^ = z'^^'e-'H^, (25) 

eq. (17) becomes Kummer's equation (Abramowitz & Stegun 1965 [13.1.1], hereafter AS65) 
or confluent hypergeometric equation (see e.g. Landau & Lifshitz 1959, Appendix d), 

zH'; + (6„ - z)H'^ - a^H^ = 0, (26) 

with a„ = 2n — 1 and 6„ = 4n — 1. The general solution of Kummer's equation is (AS65 
[13.1.11]) 

Hn{z) = AnM{an, &„, z) + 5„f/(a„, 6„, z), (27) 

where An and -B„ are arbitrary constants, and M, U are called Kummer's functions or 
confluent hypergeometric functions of the ffist and second kind, respectively (regular at the 
origin, and at infinity, respectively)^. For z = 0, M(a„,6„,0) = 1 whereas U{an,bn,z) 
diverges like z^~^" = z~^". The latter behavior implies -F„ ~ x'"" for small x, in contrast 
with the boundary condition -F„(0) = 0. Therefore Bn = and 

Fn{x) = A„(2v^)2"e-2^M(2n - 1, 4n - 1, 2y^). (28) 



^Kummer's equation describes the radial part of the Coulomb wavefunction. If a„ < 0, the solutions of 
Kummer's equation well-behaved at infinity are associated Laguerre's polynomials. 
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The constants An are fixed by imposing the boundary condition at infinity (eq. 24). Since 
M{an,bn,z) has the asymptotic behavior (AS65 [13.1.4]) 

hm M(a„,6„,z) = f " ~ ^^, 6-^"-"", (29) 

2^oo [an - 1)! 

we immediately obtain 

_ (2n-2)! 

^"-(4n-2)!^- ^^"^ 

Figure 2 shows the functions Fn{x) for n = 2-10. 

Combining the angular and radial solutions and summing over n, we finally obtain 



X- 



where 



^) = e-^^ J2 i^nX"M(2n-l,4n-l,2v^)[/iP„_i(^)-P„(/i)] (31) 

?i=2,4,... 

K -( ir/^-i 2^"(^-i)n(2n-i)! 

" ^ ^ (n-2)!!(4n-2)!(n-l)%- ^ ^ 

The function M{an,bn,z) can be evaluated numerically with standard routines, e.g. the 
Fortran program CHGM.FOR of Zhang & Jin (1996). For a; ^ 1 is better to use the full 
asymptotic expansion eq. (29), given by AS65 [13.5.1]. For a; ^ 1, the series in eq. (31) is 
dominated by the n = 2 term, and we obtain 

lim 0(x, ^) = — x^ sin^ ^, (33) 

x^o 30 

corresponding, as expected, to a uniform magnetic field with vertical field lines. Figure 3 
shows the flux function 0(x, 0) in the midplane {6 = 7r/2) given by eq. (31), and its asymptotic 
behavior given by eq. (33). At the Ohm radius the flux is reduced by a factor ~ 100 with 
respect to the asymptotic value. 

Figure 4a, b, c show the magnetic field lines in the meridional plane of the collapse region 
at different scales, computed according to eq. (31). The horizontal and vertical axis in each 
panel are the cylindrical self-similar coordinates, w = xsin^ and z = xcos^. Panel (a) 
shows the nearly uniform magnetic field inside the Ohm radius, (-R, -2)/Tohm < 1- In this 
region the series solution eq. (31) yields a good representation with the inclusion of only the 
first six terms. Panel (b) shows the magnetic field lines in the region (i?, 2;)/rohm < 10. In 
this region the series solution eq. (31) gives a good representation including the first four 
terms. Panel (c) shows the magnetic field lines in the region (-R, -2)/rohm < 100, showing 
the asymptotic convergence to the field of a split monopole at large radii. In this region the 
series solution eq. (31) gives a good representation including the first six terms. 



3. Numerical estimates 

Eq. (33) implies that the strength of the magnetic field at the center approaches a 
constant value given in dimensional form by 

$^ 

Be = ^, (34) 

307rr^hm 

where $* and rohm are given by eq. (9) and (11), respectively. Substituting these values, we 

obtain 

4G'5/2M3 ^ ^ 

Be = ^^, 35 

showing that Be scales with the electric conductivity and stellar mass as o"^M^, a result that 
may have interesting consequences for high-mass stars. 

To estimate the numerical values of the physical quantities of our model, we assume a 
given value for the uniform magnetic field Be inside the dissipation region, and we derive the 
remaining quantities in terms of Be- An estimate of the intensity of the constant magnetic 
field in the central region can be inferred from measurements of remanent magnetization in 
meteorites. Observed values range from ~ 0.1 G in achondrites to ~ 1 G in carbonaceous 
chondrites, reaching values up to ~ 10 G in chondrules. It is unclear whether this range of 
values reflects an actual variation in the strength of the magnetizing source. Chondrules have 
randomly oriented magnetizations which strongly suggest that they record magnetic fields 
that predate the accretion of the meteorites, but the measured values are the most uncertain. 
Achondrites, on the other hand, have the least complicated magnetic mineralogies, but local 
processes such as impacts may have affected their magnetization history (for reviews, see 
Stacey 1976, Levy & Sonnet 1978, Cisowski & Hood 1991 and references therein). Hereafter 
we assume Be ^ 1 G. 

From eq. (9) and (34) we obtain the Ohm radius 

^oh.n=TTr^ ^12A:V2 _^ /^ AU, (36) 



weakly dependent on the values of the parameters. 

We note that magnetic fields of strength ~ 1 G, when bent sufficiently outwards (e.g., 
by X- winds; Shu et al. 2000), can drive disk winds (e.g., Konigl & Pudritz 2000). However, 
disk winds driven from footpoints of 2 to 6 AU will have difficulty acquiring the 200 to 
300 km s~^ terminal velocities seen in high-speed jet outffows. Nevertheless, we cannot rule 
out, on the basis of these considerations, the possibility that slow disk winds co-exist with 
fast X-winds in YSOs. 



-9- 



The corresponding value for the effective resistivity rj from the expression for the Ohm 
radius (eq. 11) is 






again, weakly dependent on the numerical values of the parameters. With our fiducial values 
A^ ^ 2, M^ ^ 1 Mq, B^^ I G, we obtain rohm ~ 8.5 AU and 77 ^ 2 x 10^° cm^ s'^. This 
is larger by a few orders of magnitude than estimated values from kinetic theory of the 
microscopic ohmic resistivity in dense gas and circumstellar disks (see Sect. 6), suggesting 
that the dissipation of magnetic flux probably occurs by some anomalous diffusion process. 

The free-fall velocity at the Ohm radius is M(rohm) ~ 14 km s~^ and the gas density, 
for an accretion rate M ^ 10~^ Mq yr~^, is of the order of 10^ cm~^. This is lower than the 
often quoted value of the density at decoupling of 10^^-10^^ cm~^ (e.g., Nishi, Nakano, & 
Umebayashi 1991), because our adopted effective resistivity is larger than the conventional 
electric resistivity. 



4. Joule heating rate 

The magnetic energy annihilated per unit time and unit volume (Joule heating rate) is 

2 ^ » ^ / o^ \ 2 



^iVxBP- ^ 



4:71 IGtt" 



r sin^ 



GM. /a^y _ ,3^ 



SiT^rir^ sin^ 6 \dr J 



where we have used eq. (7) to eliminate iS(<l>). For a; <^ 1, we can approximate (f){x) with 
the asymptotic expression eq. (33), obtaining 

1 /'d(l)V sin^^ , , 

(39) 



x" 



3 sin^ e \dx 225a; ' 



Inserting this expression in eq. (38) and integrating over a sphere of radius rohm centered on 
the origin, we obtain an approximate estimate of the total energy-dissipation rate 



^^'^' - 3ooA-^ r ^- V r ^ '"' 

675Xlr]^^ * V^0/ VlO^Ocm^s- 



^ - ^.;;7T^ - 300A;M T7^ ) ( ...n \.-i ^o- (40) 



The dependence of S on the inverse fifth power of rj can be easily understood, since 
the energy dissipation rate is proportional to the resistivity times the square of the electric 
current density (oc 1]~^'^) times the volume (oc rj^) in which the current flows (eq. 38). The 
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increase of the electric current in the hmit of small resistivity (as 1]~^) suggests that the 
anomalous source of field dissipation could be associated with current- driven instabilities 
occurring when the drift speed of the charged species becomes larger than the ion's thermal 
speed, as anticipated by Norman & Heyvaerts (1985). 

With our fiducial values we obtain S ^ 3 Lq, but given the sensitive dependence of S on 
the uncertain parameter rj this number may not be very significant. What is more interesting 
is that the adopted resistivity, which is high by conventional microscopic standards (see 
Sect. 6), cannot be much lower without violating observational constraints concerning the 
total luminosity from regions of low-mass star formation (see the reviews of Evans 1999 or 
Lada & Lada 2003). For example, decreasing the value of r^ by a factor of 10, increases the 
central magnetic field to B^ ~ 10 kG, and decreases the Ohm radius to rohm ~ 0.1 AU. 
Although magnetic fields of kilogauss strength are measured on the surface of young stars, 
the energy-dissipation rate for this reduced value of rj increases to £^ ~ 3 x 10^ Lq, which 
is unrealistic for solar-mass stars. Imposing the condition that the energy-dissipation rate 
must be lower than the total accretion luminosity GM^,M/R^, we obtain a lower limit on 
the resistivity. 



Notice also that the energy-dissipation rate S scales with the stellar mass as M^, which 
is steeper than the main-sequence luminosity-mass relation, unless the resistivity increases 
significantly with increasing M*. Such a behavior runs counter to the usual notion that 
high-mass stars possess more ionizing potential than low-mass stars. Nevertheless, high 
mass stars are formed with mass accretion rates which are more than 100 times larger than 
low-mass stars (see, e.g., Osorio, Lizano & D'Alessio 1999). In these conditions of very 
high circumstellar density the ionization front will be confined to the stellar surface and the 
penetrating ionizing agent for the pseudodisks and disks will not be ultraviolet photons but 
cosmic rays and/or X-rays. 



5. Validity of the kinematic approximation 

To check the validity of the kinematic approximation (i.e. the assumption that the infall 
velocity is dominated by the gravity of the central star), we evaluate the ratio of the Lorentz 
force per unit volume in the radial direction 

F, - rJ^ X B) X Bl„ ^ j^^^^^, (42) 
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and the gravitational force per unit volume 

F, - ^. (43) 

where the density is 

P= , 2^,,| ^W' (44) 

and Q{6) yields the flattening of density contours because of pinching forces associated with 
the radial magnetic field (see Paper I). The ratio of the two forces is 

\F,\ " V^ohmy' \dxj Q{e)sm'e' ^ ^^ 

where we have used eq. (7) to eliminate S{^) and we used eqs. (12) and (14). The RHS 
is proportional to the ratio of two characteristic times, the accretion time, tacc ~ 10^ yr, 
and the crossing time of the ohmic dissipation region, tohm ~ 3 yr. The validity of the 
kinematic approximation is nevertheless ensured by the fact that the function of x and 9 on 
the right of this expression is very small. Figure Ad shows contours of the ratio IF^I/IF^I 
computed with the function Q{9) corresponding to the case i^o = 1 of Paper I. Notice that 
the force ratio is low precisely in the equatorial regions where we might expect the formation 
of a centrifugally supported disk to take place if we had included the effects of angular 
momentum in the problem. 

We also stress that the ratio of the Lorentz and gravitational forces given by eq. (45) 
depends on resistivity as t]^^. As our solution clearly shows, the non-zero resistivity of the 
gas results in a release of the field from the central protostar, as the gravitational pull of 
the central star is no longer fully available to pin the magnetic field of the central regions. 
However, the released magnetic field can be too strong for gravity to continue to win over 
the Lorentz force. If rj is sufficiently large, gravity dominates, and quasi-steady accretion 
onto the central source is possible. If rj is too small, magnetic forces overwhelm gravity, 
and the accretion region might try to explode outwards. If the coefficient of resistivity can 
reach anomalous values, the explosion outwards (when r] is small) may be coupled with re- 
implosion inwards when 7] later becomes (anomalously) large. It is interesting to speculate 
that these alternating reconnection behaviors ("flares") might correspond to FU-Orionis 
outbursts. This mechanism may be an alternate explanation to the disk thermal instability, 
possibly aided by protoplanet or protostellar companions, which has been proposed for the 
FU-Orionis phenomenon (e.g., Kawazoe & Mineshige 1993; Bell et. al 1995; Clarke & 
Syer 1996). 
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6. Comparison with other works 

According to calculations by Stepinski (1992) the electrical resistivity in the presolar 
nebula in the range of radii 3-30 AU is in the range from ~ 10^^ to ~ 10^^ cm^ s~^ if the 
grain size is ~ 1 cm, or from ~ 10^^ to ~ 10^^ cm^ s~^ if the grain size is 0.5 /xm. Clearly, the 
physical conditions in a disk are different from those of our infall model. Wardle & Ng (1999) 
have evaluated the components of the conductivity tensor for molecular gas for a variety of 
grain models, as function of the gas density. At our fiducial value of n(H2) ~ 10^ cm~^, they 
predict a ohmic resistivity rj ^ 10^^ cm^ s^^ for a standard grain-size distribution, much 
lower than our value rj ^ 10^° cm^ s^^. Therefore, the work of Stepinski (1992) and Wardle 
& Ng (1999) suggest that the effective resistivity rj of the infalling gas had better be larger 
than the microscopic electric resistivity, or severe conflicts with observational data arise. 
Given the precedent in solar physics, this result does not come totally unexpected. Unlike 
the solar case, however, the required increase in "anomalous" resistivity may not be huge, if 
the dust grains in the infalling envelopes of protostars, in contrast with those in protostellar 
disks, have not undergone too much growth. 

Desch & Mouschovias (2001) consider the dissipation of the magnetic fleld in a col- 
lapsing molecular cloud core, during the pre-pivotal stage of evolution (t < 0). Because no 
central star is present in the calculations of Desch & Mouschovias (2001), the velocity fleld 
and the density proflle clearly differ from the free-fall behavior assumed in this work. The 
equation for the evolution of the magnetic fleld has the same form as in this work, but Desch 
& Mouschovias (2001) also include in the resistivity rj also the contribution of ambipolar 
diffusion, 

^ = ^Ohm + ^AD- (46) 

The latter is given by 

IV^P _ |Bf 

IGn^-fPiPnT^ sin^ 9 4:11-/ pip„ 
where pi and p„ are the mass density of ions and neutral, respectively, and 7 is the ion- 
neutral drag coefficient. The numerical calculations of Desch & Mouschovias (2001) seem to 
suggest that the magnetic fleld approaches a steady-state conflguration towards the end of 
the run, characterized by a nearly uniform magnetic fleld of strength ^ 0.1 G over a central 
region of size ~ 20 AU. Outside this region, the fleld decreases approximately as r~^, as 
for fleld-freezing in a quasi-static isothermal envelope (Li & Shu 1996), whereas in our case 
of post- pivotal state evolution (t > 0) the fleld decreases like r~^, as for a split monopole. 
In the region of nearly uniform fleld (for density larger than ^ 10^^ cm~^), the ambipolar 
diffusion resistivity ?7ad reaches ~ 10^° cm^ s~^, and is larger by one order of magnitude 
than the ohmic resistivity (their Fig. lb). Thus Desch & Mouschovias (2001) conclude that 



^AD = ,^_., ,_,_ 2n = 7Z7— ' (47) 
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ambipolar diffusion is entirely responsible for the dissipation (or better redistribution) of the 
magnetic flux carried by the infalling gas. We note, however, that the region of uniform fleld 
in their model contains a negligible mass (~ 0.01 Mq). Once the core enters the post-pivotal 
phase of dynamical collapse, having dissipated the fleld in such a tiny fraction of the core's 
mass is of little help with respect to solving the magnetic flux problem of star formation. 

Nakano, Nishi, & Umebayashi (2002), criticized the results of Desch & Mouschovias (2001) 
noticing that they neglected the (dominant) contribution of grains to the ambipolar diffu- 
sion resistivity. Including the grain contribution, Nakano et al. (2002) found a much smaller 
value of the ambipolar diffusion resistivity than Desch & Mouschovias (2001) and concluded 
that for densities above ~ 10^^ cm~^ the fleld is dissipated by ohmic resistivity. To make 
progress, future theoretical calculations need to study the complex relationships among fleld 
dissipation, disk formation, and grain growth, as well as to include self-consistently the back 
reaction of the magnetic fleld and its momentum and energy inputs into the gas as the fleld 
is dissipated. 

An interesting question is the relative importance of ambipolar diffusion and ohmic 
dissipation in protostellar accretion flows. The ratio of the ambipolar diffusion time to the 
ohmic dissipation time is given by 

where v^ = B /[AnpnY/'^ is the Alfven speed in the neutrals, and rj„ = ijPi)^^ is the ion- 
neutral collision timescale. Notice that the ratio of timescales depends on the inverse square 
of the magnetic fleld. Thus, as the fleld is weakened, electric resistivity will ultimately 
dominate the end stages of the process. Indeed, the Ohmic term is the only one of the three 
conventionally invoked mechanisms for dissipating magnetic flelds that is linear in B in the 
induction equation (the so-called Hall term is quadratic and the ambipolar diffusion is cubic, 
see e.g. Cowling 1957). This means that Hall and ambipolar diffusion can never completely 
annihilate B; only Ohmic dissipation can do that. 

Because the induction equation (1) with only Ohmic dissipation included is linear in 
B, the superposition principle allows us to modify the approach taken in this paper and 
adapt it to other scenarios. If the process of ambipolar diffusion and ohmic dissipation 
occur in series as the gas flows inward, then, the compressed core fleld will be less than 
the ideal value due to the effect of ambipolar diffusion in the outer regions. Adopting the 
inner limit of the outer region as the outer limit of our inner calculation, everything then 
goes through as before. One would just have somewhat different coefficients for the angular 
functions in eq. (23) and, more importantly, an overall reduction of the effective flux at 
"inflnity", and thus, an increase in A^. For example, an increase of the mass-to-flux ratio by 
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ambipolar diffusion by two orders of magnitude, as found by Tassis & Mouscliovias (2005a,b), 
implies A^ ~ 200. From equation (37), tlie electric resistivity decreases by a factor of 
~ 4. Then, from equation (36), one obtains rohm ~ 1 AU. The most important reduction 
is that the reconnection luminosity decreases substantially to mean levels £ ~ O.O2L0, a 
small fraction of the bolometric luminosity, unless the reconnection events occur not steadily 
but in infrequent, powerful flares. Finally, since this paper argues that Ohmic dissipation 
operates in the innermost region, the annihilation of the elastic "core" that appears as the 
bottom-most cell in the calculations of Tassis & Mouschovias (2005a,b) may erode the base 
from which are launched outwardly propagating Shockwaves for their "spasmodic accretion 
oscillations". The point is that the inwardly advected magnetic flux past the "last zone" is 
not simply accumulated, as assumed in their calculations, but systematically destroyed. 



7. Summary and conclusions 

Assuming quasi-steady state and a spatially uniform resisivity coefficient 77, we have 
solved the problem of magnetic field dissipation during the accretion phase of star formation. 
We have adopted the velocity field determined in a previous study of the gravitational 
collapse of a magnetized cloud (Galli et al. 2006) , and we have ignored the back reaction of the 
changed magnetic topology on the fiow (using the field-freezing calculations to provide what 
we have called a "kinematic approximation"). With these assumptions, we have solved the 
problem analytically, and checked a posteriori the validity of our approximations. According 
to our solution, the magnetic field morphology changes from radial at large distances to 
asymptotically uniform approaching the origin, so that the magnetic fiux accreted by the 
central star is zero at any time. 

To determine the value of the resistivity coefficient r] we have considered the restrictions 
imposed by measurements of magnetic fields in meteorites. These constraints require an 
effective resistivity t] ^ 10^" cm^ s~^, probably several orders of magnitude larger than the 
microscopic electric resistivity of the infalling gas. Having shown that ohmic resistivity can 
dissipate enough magnetic field to solve the magnetic flux problem satisfying the available 
observational constraints, one now needs to solve the full dynamic problem of magnetic field 
dissipation and formation of a centrifugally supported protoplanetary disk in a self-consistent 
way. 
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Fig. 1. — The angular functions Gn{6) for n = 2, 4, 6, 8, and 10. 
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Fig. 2. — The radial functions Fn{x) for n = 2, 4, 6, 8, and 10. 
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Fig. 3. — The flux function in the midplane {9 = 7r/2) as function of the distance from the 
star, in nondimensional units {solid line). The dashed line shows the uniform field solution 
given by eq. (33) valid for x <^ 1. 



-20- 





0.2 0.4 0.6 0.8 1 2 4 6 8 10 

(c) (d) 

100 I I I , I I /I I /I 1/ I )■ I /I > I /I 1/ I /I 1/ I 1 I I li I I I I I I k I I I I I I I I I I 





20 40 60 80 100 0.2 0.4 0.6 0.8 1 



Fig. 4. — (a) The nearly uniform magnetic field inside the Ohm radius, (cu, 2;)/rohm < 1- 
The horizontal and vertical axis in each panel are the cylindrical self-similar coordinates, 
zu = xsinO and z = xcos9. (b) Same as (a) in the region (w, z)/rohm < 10. (c) Same as (a) 
in the region (-07, z)/rohm < 100, showing the asymptotic convergence to the field of a split 
monopole at large radii, (d) Contours of the ratio |FL|/|-Fg| (Lorentz and gravitational forces) 
for the density profile corresponding to the collapse of the Hq = 0.5 toroid (see Paper I) in the 
region (ro, z)/rohm < 1- The kinematic approximation is formally valid in the region below 
the solid curve. The values of the parameters are M^ = 0.5 Mq, M = 2 x 10^^ Mq yr"^, 
A. = 1.7. 



